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We show that the conclusions of recent paper by Flayac et al. [Phys. Rev. B 81, 045318 (2010)] 
concerning the stability of half-quantum vortices are misleading. We demonstrate the existence 
of static half-quantum vortices in exciton-polariton condensates and calculate the warping of their 
texture produced by TE-TM splitting of polariton band. 
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Half-quantum vortices are topological excitations of multicomponent condensates with combined spin-gauge 
symmetry 1 . They have been discussed^ to be the lowest energy topological excitations in exciton-polariton con- 
densates in semiconductor microcavities with unpinned linear polarization of the condensate. The half-quantum 
vortices are expected to define the Berezinskii-Kosterlitz-Thouless (BKT) transition in this system^^ and they have 
been recently observed experimentally^ 

In a recent paper— Flayac et al. studied the effects of TE-TM splitting of the exciton-polariton band on the state 
of vortices. The authors of Ref. have not found the static half- vortex solutions of quasi-equilibrium two-component 
Gross-Pitaevskii equation (GPE), and concluded that "The half- vortices are no more stationary solutions of the spinor 
Gross-Pitaevskii equations and should not affect the critical temperature of the BKT phase transition" . The goals of 
this Comment are to show that this conclusion is incorrect, to indicate the mathematical error of Ref. [B| that prevented 
the authors to establish the stationary half-vortex solutions, and to present the correct way of solving this problem. 

The order parameter of exciton-polariton condensate is the two-dimensional complex vector ip. This vector can be 
written in terms of two circular-polarization components, ip±, as 

, x + iy x - iy 

For the polariton condensate in quasi-equilibrium, the components of the order parameter ip± satisfy two coupled 
Gross-Pitaevskii equation [Eq. (10) of Ref. In the case when TE-TM splitting is present, the static vortex solutions 
can be found numerically, but one has to use the correct asymptotic behavior of ip±. In Ref. 5> it was assumed that 
the variables can be separated and, moreover, that the phases of ip± change linearly with the azimuthal angle [see 
Eq. (12) of Ref. [|. These assumptions are incorrect, and due to this incorrect substitution (12) the static half-vortex 
solutions were not found in Ref. [5|. In fact, we show below that the circular phases of the condensate order parameter 
are nonlinear functions of the azimuthal angle. 

To find the asymptotic behavior of static solutions %j>± at large distances r from the half- vortex core (i.e., in the 
elastic region), we note that in this region the polarization of the condensate becomes linear and the amplitudes of 
the components of the order parameter becomes equal, = yJn/2, where n = n(oo) is the concentration of the 
uniform condensate. Therefore, the asymptotics of tp± can be written as 



^±(r->oo,0) = (2) 

Here we denoted the circular phases as 9^fn, using the common phase angle 9 and the polarization angle rj, as defined 
in Ref. U These angles are yet unknown functions of the azimuthal angle <f>& 

To calculate the functions rj(4>) and 9(4>) we first obtain the elastic energy of the condensate in the presence of 
TE-TM splitting. Substitution of Eq. ^ into the quasi-equilibrium Hamiltonian of polariton condensate [Eqs. (3-8) 
of Ref. |j| gives 



Hcl = 2m^ J dxdy { (V??) + (W)2 + 27 



dz J \ dz J \ dz* ) \ dz* 
Here the effective mass m* and the TE-TM splitting parameter 7 are defined by 



(3) 



A = U- + -)' 7=^, (4) 
m* 2 \mi m t ) m t + mi 

where the effective mass of transverse (or the transverse-electric, TE) polaritons is m t and the effective mass of 
longitudinal (or transverse-magnetic, TM) polaritons is mi. In Eq. we also use the complex derivative 

dz 2 \dx l dy) ' 1 ' 

The equations for the angles 9 and rj are obtained by variation of H c \. In the limit r — > 00 only azimuthal derivatives 
are to be kept in these equations and we obtain 

[1-7 cos(2u)] 9" + 27 sin(2u)u'0' = 0, (6a) 
[1 + 7 cos(2w)] u" + 7 sin(2w) (l - u' 2 - 9' 2 ) = 0, (6b) 
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where 

u{4>) = ri(4>) ~ </>■ (7) 

We note that the superfluid current J that corresponds to the Hamiltonian density (J3J is given by more complicated 

expression than the usual one. In particular, for the wave function the radial and azimuthal components of the 
current are 

hn d6 . . 

J r = 7sm(2uj— — , (8a) 



hn d6 
■h = — r I 1 - 7COs(2w)] — . 



(8b) 



The equation of continuity of the current for the static solutions, divJ = 0, is reduced to the condition dJ$/d<p = 
and gives Eq. (JB^,) again. 

In what follows we denote the vortex solutions as (k, to), with the polarization and phase winding numbers k and 
to, respectively. These numbers are defined by n(2ir) — n(0) = 2nk and 0(2tt) — 0(0) = 2itm, and they can be integer or 
half-integer provided the sum k + to is an integer^ The topological charges used in Ref. [5| are l± = to =p k. Eqs. (JS^-jb) 
have simple solutions for k — 1 and an arbitrary integer to. Namely, r\ = <f> and 6 = m<f>, so that u(<p) = 0. Only these 
solutions with /_ — Z + = 2 were claimed to exist and were analyzed in Ref. [5|. The variables are indeed separated in 
this case. Moreover, the polarization field for the (1, 0) vortex ("hedgehog" ) is purely longitudinal,— and it is described 
by the radial function of the usual vortex in one-component condensate with the longitudinal mass to; . 

Apart from trivial solutions, Eqs. ([SJi,b) can be solved for any other winding numbers (A;, to). Interestingly, there 
are two qualitatively distinct type of solutions for elementary half-quantum vortices with fc,m = ±1/2. In real 
microcavities the TE-TM splitting is small and we first present the series of these solutions in the powers of 7 -C 1. 
For the (1/2, ±1/2) half-vortices the solutions are 
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sin(20) 



• sin(2</>) 



For the (—1/2, ±1/2) half-vortices the solutions are 



± 



(j) + I sin(30) + ±- sin(6</>) + 
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jj(4>) 
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<t> + J sin(30) 
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(9a) 
(9b) 

(10a) 
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The solutions to Eqs. (JSKb) can also be found numerically as shown in Fig. 1. We have used a high value of the 
TE-TM splitting parameter, 7 = 0.3, in order to illustrate the qualitative features of the behavior of polarization and 
phase angles. Analyzing the superfluid current around the vortex core we note that the streamlines are warped with 
respect to perfect circles. Physically, the warping appears due to the change of the polariton mass with polarization. 
Using Eqs. (J5^,,b) we find the streamlines to be defined by equation 



oflnr 7sin[2u(0)] 
d<j) 1 — 7 cos [2 u (</>)] 



(11) 



The warping of streamlines is shown in the inserts of Fig. 1. 

The asymptotics found above divide the solutions of the Gross-Pitaevskii equation into topologically distinct classes 
according to the values of winding numbers k and m. The half-vortex should be found by minimizing the full Gross- 
Pitaevskii Hamiltonian within a particular topological class. This solution exists since the energy is bound from below. 
This solution is static since it is a minimum of the Hamiltonian. The half-vortices can be found either by numerical 
solution of Gross-Pitaevskii equation with the required asymptotic behavior, or by other means (e.g., by variational 
method) . 

Finally, we comment on the stability of half-vortices from the other point of view. In the half-vortex core, when 
r — > 0, one of the circular components goes to zero and is singular: the order parameter behaves as re ±1( ^ so that 
the gradient is not defined at r = 0. For example, for the (1/2,1/2) half-vortex, ^_ cx re^ and ip + = const^ If 
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FIG. 1. Showing the polarization angle T](cj>) (solid blue line) and the phase angle 9{<f>) (dashed red line) of the order parameter 
of half-quantum vortices, (a) The case of (1/2, ±1/2) half- vortices, (b) The case of (—1/2, ±1/2) half-vortices. Inserts show 
the warping of streamlines (thick green lines) of the current around the half- vortex core with respect to the perfect circles (thin 
black lines). 



a half-vortex state is created by some external means and evolves according to the time-dependent Gross-Pitaevskii 
equation, this singularity in the solution will be present all the time. It is because the Gross-Pitaevskii equation is 
regular and the vortex singularities are rigid. 9 This is the reason why the estimation of life-time, given by Eq. (24) of 
Ref. d is not valid. 

To conclude, we have shown that half-quantum vortices remain stable in the presence of TE-TM splitting of 
polariton band, but their texture becomes warped — the polarization and phase angles depend nonlinearly on the 
azimuthal angle. We calculated the warping effect far away from the half-vortex core. 
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